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Abstract. Mergers of two compact objects, like two neutron stars or a neutron star and 
a black hole, are the probable progenitor of short gamma-ray bursts. These events are also 
{^JQI promising sources of gravitational waves, that are currently motivating related searches by an 

international network of gravitational wave detectors. Here we describe a search for gravitational 
waves from the in-spiral phase of two coalescing compact objects, in coincidence with short 
GRBs occurred during during LIGO's fifth science run and Virgo's first science run. The search 
includes 22 GRBs for which data from more than one of the detectors in the LIGO/ Virgo network 
were available. No statistically significant gravitational-wave candidate has been found, and a 
, parametric test shows no excess of weak gravitational-wave signals in our sample of GRBs. The 

• 90% C.L. median exclusion distance for GRBs in our sample is of 6.7 Mpc, under the hypothesis 

' of a neutron star - black hole progenitor model. 

o 
o 



1. Introduction 

Gamma Ray Bursts (GRB) are intense flashes of 7-rays traditionally divided in two classes; long 
duration bursts (duration > 2 s) which are generally associated with hypernova explosions in 
star-forming galaxies (see e.g. [H [2] and references therein) and short duration bursts (duration 
< 2 s), thought to originate primarily from the coalescence of a neutron star (NS) with another 
compact object (see e.g. [3] and references therein), either a neutron star or a black hole (BH). 
Coalescing compact objects are believed to be promising sources of gravitational waves (GW), 
currently motivating searches by a worldwide network of GW detectors. Because the time and 
sky location of a GRB is known, the search for coincident GW signals can be performed with a 
reduced threshold, increasing the sensitivity compared to untriggered searches. In the analysis 
presented here, we adopt the GRB classification scheme recently proposed by [4|, with a type-I 
GRB denoting the class of GRBs probably created by a merger (i.e. the 'short' ones), while a 
type-II GRB denotes the class of GRBs usually associated with stellar core-collapse (a 'long' 
GRB). 

This work describes a search for GW signals associated with type-I GRBs during LIGO's 
fifth science run (S5) and Virgo's first science run (VSRl). S5 took place from 4 November 2005 



to 30 September 2007 and VSRl from 18 May 2007 to 30 September 2007. S5 comprised of the 
three LIGO detectors [6j which are located at Hanford, WA and Livingston, LA. The Hanford 
site consists of two detectors, one with 4 km long interferometer arms (HI) and one having 2 km 
long arms (H2), while the detector at Livingston has arms with a length of 4 km (LI). VSRl is 
comprised of the Virgo detector (VI, [7j), located near Pisa, Italy, with 3 km long arms. 

A detection of GWs associated with a type-I GRB would provide direct evidence of the merger 
nature of such events, having large impact on GRB physics. Besides it will be possible to measure 
component masses [SIE], measure component spins [TU], constrain NS equations of state [HI [12], 
test general relativity in the strong- field regime [13] , the graviton mass [14[ [T^ and the Lorentz 
Livariance principle [TB]. Simultaneous observation of a binary's GW signal and type-I GRB will 
also enable an independent measurement of the binary luminosity distance and redshift, probing 
the relatively nearby universe's expansion, and complementing other cosmological studies (see 
e.g. |17j and references therein). 

2. Sample and data selection 

The coalescence model suggests that the time delay between the arrival of the GW and the 
subsequent electromagnetic burst, referred to as trigger time, is of order one second or less. If 
the initial spike of GRBs is created by so-called internal shocks (produced by the impact of 
faster outmoving shells on slower moving shells in the 'standard model' of GRBs (see e.g. [2] 
and references therein), the time-delay is on the millisecond scale [21 [3] or at least below 1 second 
|18j . With a more semi-analytical description of the final stage of a NS-BH merger it has been 
calculated that the majority of matter plunges on the BH within 1 second [19j and numerical 
simulations on the mass transfer suggest a timescales of milliseconds [20] or some seconds at 
maximum [21]. These considerations prompted the use of a 6s long time window (on-source) 
around the trigger time (from -5 to +1), in which we expect the GW signal to be. 

Data around this on-source segment is chosen to estimate the background and the efficiency 
of the search, split up into 324 segment^ with a length of 6 s each, not necessarily symmetrical 
around the on-source segment. To prevent a possible loud on-source signal biasing the 
background estimation, we discard 48 s on each side of the on-source segment, which reflects 
the length of the longest used template for this search; the remaining segments constitute the 
off-source segments. Considering a padding time of 72 s on both sides of the data stretch, the 
total amount of data used for the analysis is 2190 s for each GRB. To ensure the use of only 
good quality data, any 6s long segment which overlaps a data-quality veto is discarded [22] . 
Therefore the actual used number of background trials might differ from 324. 

During the whole period of S5 and VSRl, 212 GRBs in total were observed, mainly with the 
SWIFT and HETE missions, but also with the Interplanetary Network (IPN) ^65i i66i i671 168 1 169] • 
33 of them has been found to belong to the type-I GRB class. Because we require coincident data 
from at least two detectors, only 22 GRBs have been selected for the analysis (see Table [T|. If 
data from three or more detectors were available, only data from the two most sensitive detectors 
has been used for the analysis, except for GRB 070923 (see below). Two of these 22 selected 
GRBs actually had a Tqq duration larger than 2 s, but spectral features suggest their affiliation 
to the type-I GRB class. There are 9 GRBs in H1-H2 coincident time, 11 in Hl-Ll, 1 in H2-L1, 
and 1 GRB in Hl-Ll-Vl coincident time. The latter, GRB 070923, is a special case because 
it uses data from three spatially separated detectors, which all have about the same directional 
sensitivitjH to this GRB. Worth special mention is also GRB 070201, which has been already 
analysed in a high-priority search because the progenitor was possibly located in M31, a galaxy 

^ corresponding to the minimum amount of data required for a standard analysis 
^ taking into account the response of the detector in the direction of the GRB 
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Table 1. Parameters of the 22 GRBs selected for this search. Two of the GRBs in this list have 
a Tgo duration larger than 2 s, but spectral features suggest a type-I GRB affiliation. 

only ~780 kpc away. No GW signal was found and a coalescence scenario was ruled out with 
>99% confidence at that distance ^TOj, hinting to a soft-gamma repeater progenitor scenario for 
this type-I GRB in M31 [7^172]. Nevertheless, this GRB is being reanalysed in the current work 
because of the use of a lower threshold and the use of a more sensitive likelihood-ratio based 
ranking method. 

3. Search details 

The inspiral waveforms generated by coalescing NS-NS or BH-NS systems, can be reliably 
predicted using post-Newtonian perturbation theory, until the last fraction of a second prior 
to merger. These waveforms can be used in matched filtering of noisy data from GW detectors 
to identify GW candidate events. In this search, data are analysed using the standard compact 
binary search pipeline described in detail in |74] . which is indeed based on a matched-filter 
technique to filter the data against post-Newtonian approximants of the expected GW signal. 
First a discrete bank of template waveforms is set up that span a two-dimensional parameter 
space (one for each component mass) such that the maximum loss in signal to noise ratio (SNR) 
for a binary with negligible spins would be 3% [75j . While this template bank does not contain 
waveforms with spin, we verified that the search can still detect binaries with reasonable spin 
orientations and magnitudes. The template bank spans the component masses symmetrically in 
the range [1 , 40 )Mq in total mass. 

Then the data from each of the used detectors is filtered with each of the template, and a 
trigger is generated when the matched filter SNR exceeds a specific threshold. The thresholds 
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Table 2. Summary of the results for the search for GWs from each GRB. The Antenna Response 
column contains the response for each detector in the direction of the source; a value of 1 
corresponds to optimal sensitivity and a value of corresponds to no sensitivity in that direction 
(see f75] for details). An ellipsis (...) denotes that a detector's data were not used. F.A.P. is 
the false-alarm probability of the most significant on-source candidate; trials with no candidates 
above threshold are assigned a F.A.P. of 1. The last two columns show the 90% CL exclusion 
distances to the assumed merger source. 

were selected specifically for each detector combination to reflect the individual directional 
sensitivity of a detector; details can be found in [73j. We also applied two signal-based tests 
to reduce the number of background triggers. The first test {x^ test [76J) rejects triggers with 
too large differences between the expected and measured amplitude series of the signal, while 
the second test (r^ test [77]) discards triggers depending on how long the value stays above 
a certain threshold around the trigger time. Finally, coincident triggers, found to be consistent 
in their time and mass parameters from two or more detectors [78], are recorded as candidate 
events. 

4. Likelihood ratio ranking 

To obtain the search efficiency and to calculate the likelihood ratio of a candidate event, 
simulated waveform signals (software injections) are added to the off-source data and processed 
with the pipeline described in the previous section. These so-called software injections are made 
over a wide range of parameters (i.e. masses and spins of the two objects, and inclination angle 
and polarization angle of the orbital axis, describing the orientation of the binary with respect 



051114 
051210 
051211 
060121 
060313 

060427B 
06M29 
061006 
061201 
061217 
070201 
070209 

070429B 
070512 
070707 
070714 

0707 14B 
070724 
070729 
070809 

070810B 
070923 







NS-NS 
NS-BH 




















10 

distance (Mpc) 



Figure 1. The 90% C.L. exclusion distances to putative NS-NS and NS-BH progenitor systems 
for each analysed SGRB as listed in Table [TJ The median exclusion distance for a NS-NS is 
3.3 Mpc and is 6.7 Mpc for a NS-BH system. 



to the line-of-sight to earth) to assess the full efficiency. The number of recovered injections 
compared to the total number of injections determined the efficiency of the search, as a function 
of the various parameters. 

Another quantity determined by the search is the false-alarm probability; this is the 
probability that a given candidate is found by pure chance and is obtained by the relative 
number of off-source trials yielding a statistical significance equal or greater than the on-source 
measurement. With these two quantities we define a likelihood ratio of a candidate, which is 
simply the efficiency divided by the false-alarm probability. The larger this likelihood ratio, the 
more likely it is that this candidate is a real GW. 

The analysis yielded no evidence for a GW signal in coincidence with a GRB in our sample; 
see Table [2] for detailed results. With the null observations and the large number of injection 
trials we were able to constrain the distance to each GRB assuming it was caused by a compact 
binary coalescence. 

We used the approach of Feldman-Cousins [79j to compute regions in distance where GW 
events would, with a given confidence, have produced results inconsistent with our observations, 
for a given companion mass range. We quote the median exclusion distances for two mass 
ranges, representing the merger of two neutron stars (NS-NS), with a component mass in the 
range [1, 4) Mq, and for a merger of a neutron star and a black hole (NS-BH), with the black 
hole mass in the range [7, 10) Mq. The median exclusion distance for the NS-NS systems is 
3.3 Mpc while for the more massive NS-BH systems the median exclusion distance is found to be 
6.7 Mpc; more details can be found in [73]. A graphical representation of the individual excluded 
distances can be seen in Fig. [TJ 

The most significant candidate was GRB 070201 with a false-alarm probability of 6.8 %. This 
candidate was found to be just below the threshold used in the previous analysis [70], which 



showed no candidate in the on-source segment at all. Despite the low false-alarm probability, 
this candidate is consistent with the expectation from the off-source segments. 

5. Population test 

In addition to the individual detection searches, we assess the presence of a GW signal too 
weak to stand out above background separately, but which is significant when the entire sample 
of analysed GRBs are taken together. This is being done with the non-parametric Wilcoxon- 
Mann- Whitney U test [80] , to answer the question of whether both samples (significances of the 
on-source and off-source segments) are drawn from the same population. This test calculates a U 
statistic by the sums of the ranks of both samples, and compares it to the expected value if both 
samples were drawn from the same distribution. Since this statistic is approximately Gaussian 
distributed, a one-sided probability p can be calculated that either shows the two samples are 
indeed drawn from the same distribution {p ~ 50%) or contain a population of weak signals in 
the sample {p < 5%). A large value for p would hint at a problem with the analysis, since it 
implies more significant candidates in the off-source than in the on-source. 

Applying the U test, we find that the two distributions are consistent with each other; if the 
on-source and off-source significances were drawn from the same distribution, they would yield 
a U statistic greater than what we observed 53 % of the time. Therefore we find no evidence for 
an excess of weak GW signals associated with GRBs in our sample. 

6. Summary and future w^ork 

We presented a search for a merger signal in data around 22 type-I GRBs, probably created by 
a merger of two compact objects. This analysis was based on the standard matched filtering 
techniques used also for untriggered searches, but with a lower threshold. The analysis also 
used a more sensitive likelihood ratio based ranking technique to determine the significance of 
the candidate events. No GW signal in data around any type-I GRB in our sample has been 
found, and exclusion distances to each of the GRB has been computed. The median exclusion 
distance is 3.3 Mpc for a NS-NS system, and 6.7 Mpc for a NS-BH system. A non-parametric 
Wilcoxon-Mann- Whitney U test has been applied to the sample of on-source candidates, which 
was found to be consistent with the off-source background candidates. 

With the conclusion of LIGO's fifth and Virgo's first science runs, the detectors have 
undergone a series of enhancements, followed by new data taking periods, started in summer 
2009 with S6 and VSR2 runs, alternating with interruptions for further enhancements and 
commissioning activity. A fully automated online search for merger signal from type-I GRBs 
has been set up for this period, analysing the data around each trigger in the same way as it 
has been described in this paper. Regarding the analysis, several short-term improvements are 
scheduled, such as an improvement of the background estimation, a faster analysis of the data 
and the implementation of a coherent search stage; the results of this analysis will be published 
elsewhere. 

These data taking periods will last until end 2010 - mid 2011, when a shutdown of the 
detectors is planned, in order to implement advanced detector configurations |8H 182] . With 
these detectors, which will start to be operated in ~2014-2015, the sensitivity will gradually 
improve, aiming at a factor of ^ 10 improvement with respect to the initial configuration, 
increasing the observable volume of the universe by a factor of ~ 1000. By then, coincident 
detections of GWs with type-I GRBs are expected to be very likely [83], marking the beginning 
of exciting gravitational-wave physics. 
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